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PREFACE 


The  mathematical  model  investiRation  of  Tallahala  Creek  Lake  re- 
ported herein  was  authorized  by  the  Office,  Chief  of  Engineers,  U.  S. 
Army,  on  19  August  197^  at  the  request  of  the  U.  S.  Army  Engineer  Dis- 
trict, Mobile  (SAM).  An  earlier  investigation  of  Tallahala  Creek  Lake 
was  completed  in  June  1973  and  a letter  report  was  forwarded  to  SAM. 

The  model  study  reported  herein  reflects  various  structural  and  opera- 
tional changes  in  Tallaliala  Creek  Lake  not  investigated  in  the  earlier 
study. 

The  investigation  was  conducted  during  the  period  April-November 
1975  in  the  Hydraulics  Laboratory  of  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES),  under  the  direction  of  Mr.  H.  B.  Simmons, 

Chief  of  the  Hydraulics  Laboratory,  and  Mr.  J.  L.  Grace,  Jr.,  Chief  of 
Structures  Division,  and  under  the  supervision  of  Mr.  J.  P.  Bohan,  C}iief 
of  the  Spillways  and  Channels  Brandi . The  study  was  conducted  by 
Messrs.  S.  T.  Maynord  and  B.  Loftis.  This  report  was  prepared  by 
Messrs.  Maynord,  Loftis,  and  D.  G.  Fontane. 

Directors  of  WES  during  this  study  and  the  preparation  and  publica' 
tion  of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L.  Cannon,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Sl)  units  as  follows: 


Multipl 


To  Obtain 


* To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  read- 
ings, use  the  following  formula:  C = (5/9)(F  - 32) . To  obtain  Kelvin 
(K)  readings,  use:  K = (5/9)(F  - 32)  + 273.15. 
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Mathematical  Model  InvestlRatlon 

PART  I : INTRODUCTION 

Purpose 

1.  The  purpose  of  this  study  was  to  investigate  the  temperature 
structure  to  be  expected  within  the  proposed  Tallahala  Creek  Lake  and  to 
determine  a selective-withdrawal  intake  configuration  that  will  allow 
operation  to  satisfy  downstream  water-quality  objectives. 

Project  Description 

2.  The  Tallahala  Dam  will  be  located  on  the  Tallahala  Creek 
approximately  13  miles*  north  of  Laurel,  Mississippi,  and  will  have  a 
drainage  area  of  152  square  miles  (Figure  l).  An  8 ,000-ft-long  earth- 
fill  dam  will  impound  water  for  water  supply,  water-quality  management, 
and  flood  control.  The  lake  formed  by  the  dam  will  be  approximately 

8 miles  long  and  will  provide  128,360  acre-ft  of  storage  at  el  317-55** 
which  will  be  the  top  of  the  flood  control  pool.  The  storage  allocation 
will  include  675  300  acre-ft  for  flood  control  5 3658OO  acre-ft  for  water 
supply 5 and  12,it00  acre-ft  for  water-quality  control.  The  top  of  the 
conservation  pool  will  be  at  el  306.5  and  the  bottom  of  the  reservoir 
will  be  at  el  270.0.  The  surface  area  of  the  lake  at  the  top  of  the  con- 
servation pool  will  be  U,8U5  acres.  Releases  from  the  lake  will  be  made 
through  an  outlet  works  with  a 10-ft-diam  conduit  in  the  dam  and  over  a 
300-ft-long  emergency  spillway. 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (Sl)  units  is  presented  on  page  3. 

**  All  elevations  (el)  cited  herein  are  in  feet  referred  to  mean  sea 
level. 
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The  outlet  worAs  will  consist  of  an  intake  str\ictvn*e  with 
provisions  for  making  selecti ve-wl tiidrawal  ami  flood  releases.  Tlie 
selectivo-witiidrawal  releases  will  pass  tiirougii  two  10-  by  lt.75-ft  wet 
wells,  each  of  wliicii  will  have  two  3.0-  by  i».0-ft  intakes,  one  at 
el  089.0  and  one  at  el  098.0,  with  liydraulically  operated  sluice  gates. 

The  two  intakes  at  el  098.0  will  discharge  into  1,he  aides  of  the  wet 
wells,  and  the  two  at  el  289. 0 will  discharge  into  the  front  faces.  The 
flows  fi-om  each  of  ttie  wet  wells  will  be  conti'olled  by  2.0-  by  0.5-ft 
sluice  gates  at  ol  2T0.0.  The  flows  from  the  wet  wells  will  discharge 
into  the  flood  control  conduit  Just  downstrefim  of  the  service  gates. 

Fneh  wet  well  will  pass  85  cfs,  and  the  combined  capacity  for  both  wet 
wells  will  be  ITO  cfs. 

ApprOvUch 

!».  The  study  was  accomplished  with  the  use  of  a numerical  siimila- 
tion  model.  The  approach  involved  the  selection  of  several  study  years 
and  simulation  of  lake  operation  for  each  of  these  years.  St.udy  years 
selected  had  combinations  of  stre:\mflow  quantities  and  air  temperatures 
that  could  create  extreme  conditions  of  thermal  sti-ati f ication . The 
data  required  for  the  simulations  were  lake  inflows  and  outflows,  inflow 
stresun  temperatures,  meteorological  data  for  each  of  the  study  years, 
geometry  of  the  lake,  and  geometry  of  the  intake  structure. 

5.  The  heat  transfer  into  and  out  of  the  lake  was  evaluated  and 
the  heat  was  distributed  wit.hln  the  lake.  A heat  budget  was  maintained 
throughout  the  simulation  period.  An  objective  temperature  was  speci- 
fied for  each  simulation  day.  und  an  operating  scheme  was  determined. 

The  operation  scheme  for  any  day  was  the  combination  of  open  ports  that 
minimized  the  difference  between  the  objective  downstroiun  temperature 
and  the  predicted  release  temperature.  The  output  from  the  simulation 
included  a comparison  between  objective  and  release  temperatxire  in  graph- 
ical form  throughout  the  simulation  period  as  well  as  tabular  summaries 

for  each  day  and  plotted  profiles  of  temperature  witliln  the  lake  at  j 

specified  times  of  the  year.  The  tuunerical  simulation  model  was  also 
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used  to  simuliite  the  dissolved  oxygen  (D.O.)  structure  of  Tnllahala 
Creek  Luke.  The  purpose  of  these  simulations  was  to  evaluate  the  ef- 
fects of  reservoir  operation  upon  the  D.O.  regime  of  the  lake.  A de- 
scription of  the  simulations  tuid  the  results  of  the  simulations  are  pre- 
sented in  Appendi.x  A of  this  report. 
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PART  II;  MODEL  DESCRIPTIOH 


6.  The  downstream  release  characteristics  and  the  internal  temper- 
ature and  D.O.  structure  for  Tallahala  Creek  Lake  were  predicted  using 

a numerical  simulation  model.  The  model  (WESTEX)  used  in  conjunction 
with  this  investigation  was  developed  at  the  t'.  S.  Army  Engineer  Water- 
ways Experiment  Station  (WES)  based  on  the  results  of  Clay  and  Fruh,^ 

2 3 L 

Edinger  and  Geyer,  Dake  and  Tarleman,  and  Bohan  and  Grace.  The  D.O. 

portion  of  the  WESTEX  model  is  described  in  Appendix  A. 

7.  The  WESTEX  model  provides  a procedure  for  examining  the  bal- 
ance of  thermal  energy  imposed  on  an  impoundment.  This  energy  balance 
and  lake  hydrodynamic  phenomena  are  used  to  map  vertical  profiles  of 
temperature  in  the  time  domain.  The  model  includes  computational  meth- 
ods for  simulating  heat  transfer  at  the  air-water  interface,  heat  advec- 
tion  due  to  inflow  and  outflow,  and  the  internal  dispersion  of  thermal 
energy.  The  model  is  conceptually  based  on  the  division  of  the  impound- 
ment into  discrete  horizontal  layers.  Fundamental  assumptions  include 
the  following: 

a.  Isotherms  are  parallel  to  the  water  surface  both  laterally 
and  longitudinally. 

b.  The  water  in  each  discrete  layer  is  isotropic  and  physi- 
cally homogeneous. 

£.  Internal  advection  (between  layers  in  the  lake)  and  heat 
transfer  occur  only  in  the  vertical  direction. 

External  advection  (inflow  into  and  outflow  from  the  lake) 
occurs  as  a uniform  horizontal  distribution  within  each 
layer. 

e.  Internal  dispersion  of  thermal  energy  is  accomplished  by  a 
diffusion  mechanism  which  combines  the  effects  of  molecu- 
lar diffusion,  turbulent  diffusion,  and  thermal  convection. 

8.  The  surface  heat  exchange,  internal  mixing,  inflow,  and  out- 
flow processes  are  simulated  separately  and  their  effects  are  introduced 
sequentially  at  daily  intervals. 

9.  The  WESTEX  model  employs  an  approach  to  the  evaluation  of 
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surface  heat  transfer  that  was  developed  by  Edinger  and  Geyer.  This 
method  formulates  equilibrium  temperatures  and  coefficients  of  surface 
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f'imalat.ion  rtf  t.ho  proee;-!;-,  of  inflow  ti  i apl  aot'r.  \ipw!ird  a vttlame  t-tpial  to 
till'  total  inflc'w  qajint.ity.  This  apwarii  ti i !iplaeement.  is  refleeleti  in  the 
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model  by  an  incronao  in  ttie  water  luirt'ace.  A eorreapoiul i tiK  derreaae  in 
water  r.urt'ace  occurs  as  a result  of  t.lie  outflow  process. 

iH.  The  int.ernal  dispersion  process  is  represent('d  by  an  internal 
mixiny:  scheme  hasiai  on  a simple  diffusion  analor^v.  Internal  mixing 
transfers  heat  and  other  wator-ipial  i t.y  conc.t  i t uent.s  between  adjacent 
layers.  The  magnitude  of  the  ti'anc.fer  between  two  layers  is  a pei'cent- 
ay.e  of  the  total  tiainsfer  required  to  completi'ly  mix  the  two  layers. 

This  percenta^’:e  is  a mixirut  coefficient  that  is  liefiin'd  for  every  layer. 
Data  input  includes  values  of  the  mixing  coefficient  at  the  top  and  at 
the  bottom  of  t,ht'  lake.  An  exponential  fit  betwt'cn  the  t.wo  extreme 
values  is  used  to  determine  the  appropriate  coefficient,  at  each  layer. 

15.  The  outflow  component  of  the  model  incorporates  the  selective- 

)j 

withdrawal  techniques  developed  I'.v  Itohan  and  ilrace.  Transcendental 
equations  defining  t lie  r.ero  velocity  limits  of  (he  withdrawal  zone  are 
solved  witli  a half-interval  search  method.  With  knowledge  of  the  wit.li- 
drawal  limits,  th('  velocity  profile  due  to  outflow  can  be  determined. 

The  flow  from  each  layer  is  then  the  product  of  the  velocity  in  the 
layer,  the  widtli  of  tlie  layer,  and  the  thickness  of  the  layer.  A flow- 
weighted  average  is  applied  to  water-qual i t.y  profiles  to  determine  the 
value  of  the  release  content  of  each  pariuneter  for  each  time  step. 

16.  The  lake  regulation  algorithms  liave  been  developed  to  realis- 
tically simulate  the  field  operation  of  a selective-withdrawal  system. 

The  selective-withdrawal  system  is  assumed  to  be  configurini  with  an 
arbitrary  luunber  of  selective-withdrawal  intakes  located  in  <’ach  of  (wo 
wet  wells  with  a separate  floodgate.  Maximiuii  flows  and  minimum  flows 
from  each  intake  and  from  t.he  floodgate  must  be  specified.  Also,  the 
maximum  flow  for  the  sel  ect.i  ve-wi  thdrawal  syst.em  is  specified.  The 
algorithms  at.tianpt.  to  numerically  withdraw  water  at.  or  near  the  objec- 
tive temperature.  Withdrawal  will  be  from  either  one  intake  l<'vel,  two 
adjacent  intake  levels,  and/or  the  flood  control  intake  depending  upon 
the  objective  temperat.ure,  the  tempi'rat.ure  (irofile,  tlu'  intake  capaci- 
ties, and  the  .'uiiount  of  flow  t.o  be  released. 
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PART  III:  DFRTELOPMENT  AND  ACQUISITION  OF  DATA 


r 


Selection  of  Study  Years 


IT.  For  the  selection  of  study  years,  statistical  analyses  of 
mean  monthly  strefimflow  and  mean  monthly  di-y  bulb  temperature  were  con- 
ducted for  the  jieriod  of  record  19^8-1968  (Plate  l).  Study  years  were 
limited  to  this  period  due  to  lack  of  adequate  meteorological  data  prior 
to  191*8  and  lack  of  stretvmflow  records  after  I968.  Only  records  from 
March  through  October  were  considered  in  the  selection  of  study  years. 
Experience  has  shown  that  this  is  the  period  in  which  density  strati- 
fication in  the  lake  is  affected  by  hydrology  and  meteorology.  Emphasis 
was  given  to  the  characteristics  of  the  sjiring  months  due  to  the  partic- 
ular importance  of  these  months  in  fish  reproductive  cycles. 

18.  Combinations  of  above  average,  average,  and  below  average 
t\,vdrologic  and  meteorologic  conditions  were  considered  in  the  selection 
of  study  years.  The  five  years  discussed  below  were  selected  for  the 
analysis  of  temperature  at  Tallahala  Creek  Lake. 

a.  I95I*  - Below  normal  runoff  occ\irred  throughout  the  year. 
Air  temperatures  were  above  normal  in  February,  April, 
and  June  tdirougii  September.  This  condition  would  tend  to 
allow  stratification  to  form  early  in  the  year  and  remain 
well  into  the  fall . 

1957  - Very  low  flows  occurred  in  the  beginning  of  tiie 
year  and  were  accompanied  by  above  normal  ai?-  tempera- 
tures. The  summer  period  had  below  normal  runoff  with 
average  air  temperatures;  above  normal  runoff  oceurred 
from  September  through  December.  Stratification  would 
tend  to  form  early  in  the  year  and  then  liecay  early  due 
to  tile  iiigh  flows  in  the  fall . 

£.  1958  - The  runoff  was  normal  tlirougliout  tlic  year,  and  itie 

air  temperatures  were  below  normal  from  January  througii 
March  and  normal  from  April  througii  November. 

1961  - Well  above  normal  runoff  occurred  from  February 
througii  April  and  in  December  with  nearly  normal  flows 
for  the  remainder  of  tlie  year.  Air  t.emperatures  were 
below  normal  in  January  and  from  April  througii  October. 
Tiie  iiigii  spring  flows  would  iniiibit  tiie  formation  of 
stratification,  and  tiie  low  air  temperatures  during  tiie 
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summer  would  result  in  lower  tVian  normal  water  tempera- 
tures near  the  surt’ace. 

£.  1961*  - March,  April,  and  December  had  above  normal  runoff. 

The  remainder  of  the  year  had  nearly  noi’mal  flows.  The 
monthly  air  temperatures  were  nearly  normal  with  the  ex- 
ception, of  April  (above  normal  ) and  Februar’y  and  October 
(below  normal).  The  high  flows  in  March  and  April  would 
tend  to  inhibit  the  formation  of  strong  stratification, 
wliile  conditions  throughout  the  remainder  of  the  year 
were  nearly  normal. 

Data  Requirements 

Meteorology 

19*  Meteorological  data  from  the  Meridian,  Mississippi,  Weather 
Station  were  used  for  this  study.  The  weather  station  is  locat.ed 
1*0  miles  north  of  the  diimsite.  The  required  data  consist  of  dry  bulb 
temperature,  dew  point  temperature,  wind  speed,  and  cloud  cover.  These 
data  wore  obtained  from  the  National  Climatic  Center  in  Asheville,  North 
Carolina.  Eight,  observations  were  furnistied  for  each  day.  Daily  aver- 
51(^0  values  were  computetl  and  used  to  determine  equilibrium  temperatures, 
surface  heat  e.xchange  coefficients,  and  daily  average  net  solar  radi- 
ation quantities  for  the  period  of  record. 

Hydrology 

20.  Mean  daily  inflow  and  outflow  qtiantities  are  shown  in  Elate  2. 
Hydrologic  routings  were  conducted  by  the  U.  S.  Army  F:ngineer  District, 
Mobile  (SAM),  to  determine  these  flows. 

Lake  geometry 

21.  The  area-volume  curve  is  shown  in  Elate  1.  Tliis  curve  and 
other  data  describitig  the  location  and  design  of  tlie  intake  structure 
were  furnished  by  SAM. 

Stream  temperature 

22.  Stream  temperature  recoi'ds  for  the  five  stiuiy  years  were  not 
available.  Some  stream  temperatures  were  measured  between  February  and 
August  1965  at  iiaurel  , Mississippi,  approximately  10  miles  below  the  d!un- 
site.  These  temperatvires  were  reported  in  t.he  Eascagoula  River  Compre- 
hensive Basin  Study.  These  data  wei'e  used  in  t.he  development  of  a 
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rogression  equation  relating  equilibrium  temperature,  streamfiow,  aiui 
observed  stream  temperature.  The  following  regression  model  was  used. 


a + 


'i^t " 


‘'3''t-l  " 


(M 


where 


and 


0 = mean  daily  strejim  temperat.ure , °F 
t = Julian  diiy 

Q = mean  daily  streamfiow,  efs 
E = mean  daily  equilibriiun  temperature,  ”F 
ct  , 3 are  regressioti  coefficients  as  follows: 

rt  = 19.'»T5 


3j  = -0.0020 

3^  = 0.13595 
3^  = 0.1231* 

3)^  = 0,1*095 

Tins  equation  has  a correlation  coefficient  of  0.905  and  a standard  error 
of  2.21°F.  The  mean  daily  lake  inflow  temperatures  (Plate  h)  for  the 
five  stiuly  years  were  predicted  using  this  equation, 

Ob.lective  tempei'ature 


23.  A least-squares  analysis  was  used  to  fit  a harmonic  curve  to 
the  predicted  stream  temperatures  for  the  five  study  years.  The  curve 
represents  the  .average  natural  stream  temperature  variation  to  be  ex- 
pected during  !i  year.  The  following  regression  model  was  used: 


0^  = A sin  (Pt  + C)  + D (5) 

The  coefficient  B is  a unit  conversion  from  diiys  to  radians.  The  co- 
efficients A , C , and  D were  determined  by  solution  of  Equation  5 
with  the  Newton-Raphson  technique  and  were  computed  to  be  the  following: 
A = -13.831* 

B = 1.721  X 10“^ 

C = 1.333 
D = 65.1*1* 


ll* 


> 
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Equation  5 was  used  to  defitjo  the  dowustro.'un  temperature  ol'Jeetive. 

2h.  The  entire  record  of  predicted  downstroun  temperatures  was 
scanned  foi-  ttie  maximum  stream  temperature  to  be  expected  for  each  day 
of  ttie  year.  These  365  maximum  t.emt'eratures  were  ttien  fit  to  t)ie  name 
rottression  model  as  indicated  in  Equation  5.  A similar  sine  curve  was 
determined  for  tlie  mininuuti  temperatures  to  be  expected  cacti  day  of  the 
year  over  ttie  period  of  record.  The  coefficients  for  these  curves  are 
as  follows; 


Coefficient 

A 

B 

C 

D 


Maximum 

-11 .91* 

1 .721  X 10~^ 

1.325 

70.63 


Mini mum 

-15.66 
1.721  X 10 
i . 325 
60.25 


These  ciu'vos  of  maximum  and  minimum  predicted  dowTintremii  temperatures 
(■tive  an  indication  of  tlie  variation  of  natural  stremn  temperatures  from 
tile  single  harmonic  curve  that  is  used  as  an  objective  temperature. 


Model  Calibration 


25.  As  has  been  discussed  previously,  the  WEGTEX  model  requires 
the  determination  of  coefficients  of  surface  heat  exchange  distribution 
and  internal  mixint';.  Eor  Tallahala  Creek  Lake  these  coeff icient.s  were 
detennined  by  conductitiK  simulations  with  Tallahala  hydrologic  and  mete- 
orologic  data.  Coefficients  were  adjusted  and  simulation  was  repeated 
until  the  predicted  temperature  profiles  corresponded  in  shape  and  range 
to  those  observed  in  nearby  Okatibee  I.ake.  Okatibec  hake  is  located 
Just  north  of  Meridian,  Mississippi.  It  is  similar  in  size,  dept.li,  and 
flow  mjignitudes  to  Tallahala  Creek  Lake.  I’rofiles  of  temperature  and 
D.O.  in  Okatibee  I.ake  were  obtained  from  GAM  (Plate  5).  Isotherms  pre- 
dicted for  Tallahala  Creek  Lake  (Plate  8)  compare  favorably  wi  t.ti  Okat  ibee 
data.  The  following  coefficients  were  determined  from  the  analysis: 

B = 0.6 

A = 0.2 
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= 0.5 


= 0.3 

where 

B = percentage  of  incoming  shortwave  radiation  absorbed  in  the 
surface  layer 

\ = absorption  coefficient,  ft~^ 

= mixing  coefficient  at  surface 
= mixing  coefficient  at  bottom 

Since  data  did  not  exist  to  accurately  determine  the  amount  of  entrain- 
ment induced  by  Tallahala  Creek  Lake  inflows,  simulations  were  conducted 
assuming  no  entrairunent  and  entrainment  of  a volume  equal  to  the  volume 
of  inflow.  The  simulati''ns  indicated  that  the  effect  of  entrainment  on 
the  predicted  thermal  profiles  and  predicted  release  temperatures  were 
negligible,  and  no  entrainment  by  inflow  was  assiuned  for  all  subsequent 
simulations. 

26.  As  mentioned  previously,  two  inlets  in  the  proposed  intake 

structure  will  discharge  into  the  sides  of  the  wet  wells.  A previous 

7 

physical  hydraulic  model  study  on  the  outlet  works  for  Mew  Hope  Res- 
ervoir (now  designated  B.  EVerett  Jordan  Reservoir)  conducted  at  WEC 
had  indicated  that  side  inlets  could  have  different  selective-withdrawal 
characteristics  compared  with  inlets  on  the  front  face  of  the  intake 
structure.  For  the  same  discharge  conditions,  the  inlets  located  on  the 
upstream  face  of  the  New  Hope  intake  structure  were  approximately  twice 
as  effective  in  withdrawing  water  from  above  the  density  interface  as 
were  those  located  on  the  sides.  For  the  Tallahala  intake  structure  it 
was  thought  that  the  two  proposed  side  inlets  would  not  have  different 
selective-withdrawal  characteristics  from  front  facing  inlets.  The  New 
Hope  intake  structure  was  recessed  back  into  the  earth-fill  dam  such 
that  the  fill  caused  shallow  depths  at  the  side  inlets.  Also,  the  local 
topography  at  New  Hope  caused  restricted  access  of  flow  to  the  side  in- 
lets in  the  structure.  The  combination  of  these  two  effects  was  be- 
lieved to  be  responsible  ^or  the  different  selective-withdrawal  charac- 
teristics of  the  side  inlets.  The  intake  structure  at  Tallahala  Creek 
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Lake  extends  t’ai'tlier  into  the  reservoir,  and  site  conditions  are  not  re- 
strictive compared  wi  t}i  the  New  Hope  structure.  Tn  addition,  i.herc  is  a 
si^^ni ficant  difference  in  the  discharges  to  he  released  through  the 
Tallaiiala  and  New  Hope  structures.  Discharges  tested  in  the  New  Hope 
stuiiy  ranged  from  300  to  2700  cfs  wliereas  the  selective-withdrawal 
capacity  of  tlie  Tallahala  stinicture  is  175  cfs.  Therefore,  for  this 
stvuiy,  the  selective-witiuirawal  technique  in  the  WIOSTKX  model  was  as- 
sumed to  apply  equally  to  side  or  front,  facing  inlets. 
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PAJ^T  IV:  SIMIIIJVTIONS 


27.  Four  configurations  were  used  in  the  analysis  of  the  location 
and  operation  of  the  selective-withdrawal  intakes.  Three  of  the  config- 
urations have  selective-withdrawal  capabilities.  The  fourth  configura- 
tion allows  floodgate  releases  only. 

28.  The  capacity  of  the  selective-withdrawal  system  was  estab- 
lished such  that  the  system  would  be  large  enough  to  pass  most  flows. 
Outflow  requirements  larger  than  the  selective-withdi’awal  capacity  must 
be  satisfied  by  operation  of  the  floodgates.  Plate  6 shows  an  outflow- 
exceedance  curve  that  was  determined  from  a statistical  analysis  of 
daily  outflow  requirements  from  April  to  September  for  all  of  tiie  study 
years.  The  exceedance  curve  shows  that  for  possible  selective-withdrawal 
capacities  above  about  130  cfs,  even  a large  increase  in  capacity  will 
not  significantly  increase  the  percentage  of  flows  which  can  be  con- 
trolled selectively.  For  Tallahala  simulations,  the  selective-withdrawal 
system  capacity  was  established  by  SAM  to  be  170  cfs.  This  capacity  al- 
lows selective  control  of  85  percent  of  the  flows  for  the  study  years. 
Thus,  operation  of  the  floodgate  is  required  for  15  percent  of  the  re- 
lease flows  from  April  to  September  foi’  the  study  years. 

29.  The  selective-withdrawal  system  was  simulated  as  being  con- 
figured with  two  wet  wells  each  containing  two  intakes.  Wliile  it  is 
hydraulically  possible  to  release  flow  through  two  Intakes  of  the  stune 
wet  well,  this  practice  is  not  recommended  for  it  may  result  in  unequal 
flow  distribution  from  the  two  intakes  and  ineffective  blending  as  well 
as  induce  unstable  hydraulic  flow  conditions.  If  blending  of  flows  is 
required  to  meet  a temperat\ire  objective,  the  model  will  use  one  intake 
location  from  each  wet  well. 

30.  The  first  intake  conf iguratioti  (tyqie  l)  considered  consisted 
of  two  wet  wells  and  a floodgate.  Eacti  wet.  well  contained  selective- 
withdrawal  intakes  with  center  lines  at  el  298.0  and  289.0.  Each  intake 
had  a capacity  of  85  cfs.  The  flood-release  conduit  was  a 10-ft-diam 
circular  conduit  with  invert  at  el  272.0.  Opei’ation  of  this  system  al- 
lowed simultaneous  releases  from  the  floodgate  and  selective-withdrawal 
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intakes.  of  floodgate  releases  witli  selective-withdrawal  re- 
leases was  used  to  achieve  a downstream  temperature  ob.lective  during  i 

times  when  cold  water  was  needed.  Tliis  blending  of  floodgate  releases 
might  be  accomplished  through  the  use  of  a low-level  outlet  witliin  the 
floodgate. 

31.  The  second  configuration  (type  II ) considered  was  proposed 
by  SAM  (see  paraj^^raph  3).  It  consisted  of  two  wet  wells  and  a floodgate 
Kach  wet  well  had  selective-withdrawal  inttxkes  with  center  lines  at 
el  298.0  and  280. 0.  Each  intake  had  a capacity  of  85  cfs  and  a minimum 
controllable  flow  of  2 cfs.  The  flood-release  conduit  was  a lO-ft-disim 
circular  conduit  with  invert  at  el  272.0.  Operation  of  this  system  per- 
mitted floodgate  release  only  when  the  required  ovitflow  was  greater  than 
the  capacity  of  the  selective-withdrawal  system. 

32.  The  third  configuration  (type  III)  considered  consisted  of 
two  wet  wells  and  a floodgate.  One  wet  well  contained  selective- 
withdrawal  intakes  witli  center  lines  at  el  298.0  and  286.0  and  the  otlier 
wet  well  contained  int;ikes  witti  centei-  lines  at  el  301*. 0 and  202.0. 

Each  intake  had  a capacity  of  85  cfs  and  a mininuun  cont.rollable  flow  of 
2 cfs.  Tlie  flood-release  conduit  was  a 10-ft-diajii  circular  conduit  witli 
invert  at  el  272.0.  Operation  of  this  system  permitted  blending  of 
selective-withdrawal  and  floodgate  releases.  The  primary  advantage  of 
type  III  over  type  II  was  to  allow  withdrawal  of  warm  water  from  a 
higher  pool  elevation  during  the  period  in  which  warm  water  was  required 
for  the  downstreiun  temperature  object, ive.  Also,  as  with  type  I,  more 
control  was  allowed  for  release  of  cold  water  due  to  the  availability 
of  blending  floodgate  releases  with  selective-withdrawal  releases. 

33.  The  fourth  configuration  (type  IV)  considered  consisted  of  a 
floodgate  only.  The  flood-release  conduit  is  a 10-ft-dijim  circular  con- 
duit with  invert  at  el  272.0.  Operation  of  this  system  required  flood- 
gate releases  for  all  outflows.  The  purpose  of  considering  all  releases 
to  be  through  tlie  floodgate  was  to  assess  the  effect  on  dowiistream  re- 
lease temperatures  if  no  selecti ve-withdi'awal  facilities  are  provided. 

It  is  assumed  that  thin  approximates  a worst -case  operation  condition 


since  the  presence  of  selective-withdrawal  facilities  allows  control  for 
warmer  releases. 


3U.  Simulations  were  conducted  for  intake  configurations  types 
I-IV  for  all  five  study  years.  Plate  7 shows  computed  release  temper- 
atures and  the  required  daily  intake  hydrographs  for  all  simulations. 

The  predicted  structure  of  temperature  within  the  lake  during  each  of 
the  five  study  years  is  shown  in  Plate  8.  The  isotherm  plots  reflect 
results  from  simulations  with  a type  II  intake  configuration.  Simula- 
tions with  types  I and  III  intake  configurations  yielded  similar  iso- 
therm plots. 

35.  Additional  simulations  were  conducted  with  a simple  operation 
plan.  The  description  and  results  of  these  simulations  are  given  in 
Appendix  B.  These  simulations  indicated  that  a simple  operation  plan 
could  meet  temperature  objectives. 
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PART  V:  DISCUSSION 


36.  An  attempt  was  made  in  this  study  to  satisfy  an  objective 
temperature  band  based  on  the  maximum  and  minimum  natural  stream  temper- 
atures for  the  five  study  years.  In  general,  the  release  temperatures 
fell  within  the  objective  band  for  intake  configurations  with  selective- 
withdrawal  capability  (types  I-III ) . No  appreciable  differences  exist 
between  the  results  for  the  five  study  years. 

37.  Results  of  simulations  of  the  intake  configuration  without 
selective-withdrawal  capability  (type  IV ) yielded  computed  release  tem- 
peratures that  were  cooler  than  natural  stream  temperature  during  the 
spring  months.  This  is  the  period  during  which  cooler  than  natural 
stream  temperatures  could  adversely  affect  fish  spawning. 

38.  For  results  of  the  simulations  with  selective-withdrawal 
capability  (types  I-III ),  computed  release  temperatures  from  mid-August 
through  mid-October  are  above  the  objective  temperature  band.  Analysis 
of  increased  selective-withdrawal  capacity  (type  II,  300  cfs)  and  low- 
level  withdrawal  capabilities  (represented  by  floodgate  blending  in  the 
types  I and  III  configurations)  yielded  no  significant  improvement  in 
the  results.  The  failure  of  the  project  to  meet  the  objective  band  in 
certain  instances  was  a result  of  either  insufficient  available  storage 
of  desired  temperature  water  or  floodgate  operations  necessitated  by 
very  large  outflows. 

39-  The  results  of  the  simulations  show  minor  differences  between 
the  intake  configurations  containing  selective-withdrawal  capabilities 
(types  I-III ).  It  appears  that  the  type  II  configuration  would  be  ad- 
vantageous since  it  would  be  simpler  to  operate.  It  is  recommended  that 
the  type  II  configuration  be  used  and  operated  in  accordance  with  the 
simple  operation  plan  (Appendix  B).  It  is  expected  that  the  proposed 
outlet  works  will  provide  sufficient  reaeration  to  achieve  acceptable 
levels  of  D.O.  (5  mg/Jl  minimum)  in  the  releases  regardless  of  the  D.O. 
content  of  the  flow  entering  the  outlet  works  and  plan  of  operation 

I (Appendix  A). 
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APPENDIX  A:  DISSOLVED  OXYGEN 


1.  The  dissolved  oxygen  (D.O.)  content  of  the  reservoir  can  be 
defined  as  the  sources  minus  the  sinks.  The  sources  of  D.O.  in  impound- 
ments include  photosynthesis  in  the  euphotic  zone,  atmospheric  reaera- 
tion, and  reservoir  inflows.  D.O.  sinks  include  respiration  of  plants 
and  animals,  oxygen  demand  of  reservoir  inflows,  benthic  and  detritus 
oxygen  demands,  and  D.O.  in  reservoir  outflows. 

2.  Observations  of  actual  D.O.  profiles  in  various  lakes  by 

8*  9 10 
Bella,  Carroll  and  Fruh,  and  Fontane  and  Bohan  indicate  that  a por- 
tion of  the  lake  below  the  surface  can  experience  temperature-dependent 
saturated  D.O.  conditions.  For  depths  below  this  saturation  zone,  the 
net  effect  of  all  D.O.  sources  and  sinks  can  be  represented  mathemati- 
cally by  a total  D.O.  depletion  term. 

3.  The  WESTEX  model  contains  a simple  method  for  routing  D.O. 
based  on  the  work  of  Bella, ^ Carroll  and  Fruh,^  and  Fontane  and  Bohan. 
D.O.  is  routed  in  a manner  similar  to  that  used  for  temperature.  The 
processes  simulated  in  the  model  are  advection,  surface  saturation, 
internal  dispersion,  and  oxygen  depletion.  The  D.O.  routing  portion  of 
the  WESTEX  model  is  used  to  evaluate  the  relative  effects  of  structural 
design  and  project  operation  on  the  D.O.  budget  of  the  lake.  The  re- 
sults should  not  be  interpreted  as  predictions  of  absolute  day-to-day 
D.O.  concentrations.  As  indicated  by  Bella, ^ the  resultant  D.O.  pre- 
dictions do  not  account  for  short-term  D.O.  variations  (for  example,  due 
to  an  algal  bloom),  but  rather  reflect  the  progressive  D.O.  changes  that 
occur  with  depth  over  the  entire  stratification  period. 

Advection  and  Internal  Dispersion 

The  D.O.  content  of  the  inflow  and  the  outflow  is  evaluated  and 
used  to  adjust  profiles  within  the  lake.  An  internal  mixing  mechanism. 


* Raised  niunerals  refer  to  similarly  numbered  items  in  the  References 
at  the  end  of  the  main  text. 
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based  on  an  analogy  of  the  thermal  diffusion  process,  is  applied  to  the 
D.O.  profile.  This  internal  mixing  mechanism  simulates  the  movement  of 
D.O.  from  the  upper  to  tlie  lower  layers  of  the  lake. 


Oxygen  Depletion 

5.  Total  oxygen  depletion  rates  are  determined  by  plotting  D.O. 
versus  time  for  a given  elevation.  This  depletion  rate  is  a function  of 
all  the  processes,  hydrodynamic  and  oxygen-consuming,  that  occur  in  the 
reservoir.  The  WESTEX  model  accounts  for  the  hydrodynamic  processes; 
therefore,  the  val\ie  of  the  depletion  rate  used  in  the  WESTEX  model  is 
selected  to  represent  only  the  oxygen-consuming  processes.  The  de- 
pletion rate  used  in  the  model  is  selected  such  that  the  D.O.  profiles 
produced  by  the  model  will  yield  total  depletion  rates  which  are 
reasonable  in  terms  of  observed  data  on  similar  reservoirs. 

6.  Within  the  WESTEX  model,  the  depletion  rate  for  D.O.,  a con- 
stant value  throughout  the  entire  simulation  period,  is  applied  each 
daily  time  step  to  every  layer  below  the  surface  saturation  zone. 

Surface  Saturation 

7.  As  indicated  previously,  observed  D.O.  profile  data  liave  shown 
the  existence  of  a surface  saturation  zone.  Many  factors  influence  the 
characteristics  of  this  zone.  Considerable  analysis  of  many  lake  pro- 
files is  needed  for  a general  determination  of  saturation  zone  thickness 
and  percentage  of  D.O.  saturation  within  this  zone.  The  values  used  for 
a given  project  study  should  be  based  on  observed  data  from  similar 
lakes. 

Dissolved  Oxygen  Data 

3.  Data  on  the  D.O.  content  of  flow  in  Tallaliala  Creek  were  avail- 
able from  the  preimpoundment  study^^  of  the  project.  Based  on  these 
data,  two  regression  equations  were  developed  to  predict  tlie  D.O. 
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content  of  the  inflow  to  Tallahala  Creek  Lake: 


and 


D.O..  = -2.88  + l.lU  D.O.  ^ for  T > T°C 

inflow  sat 


D.O. 


Inflow 


0.90  X D.O. 


sat 


for  T < 7°C 


where 

D.O. . = inflow  D.O.  content,  mg/H 

inflow 

T = stream  temperature,  °C 

D.O. sat  ~ temperature-dependent  saturated  D.O.  content,  mg/?. 

The  saturated  D.O.  content  was  computed  by 

^■°‘sat  0.0677  + O.OO208T 

Daily  values  of  the  computed  D.O.  content  of  tlie  inflow  are  shown  in 
Plate  A1  for  each  of  the  five  study  years.  The  net  D.O,  contributed  by 
inflow  is  the  inflowing  D.O.  content  decreased  by  an  amount  whicli  repre- 
sents D.O.  depletion  due  to  the  travel  time  within  the  lake  required  for 
the  inflow  current  to  reach  the  dam.  The  travel  time  for  an  inflow  cur- 
rent to  reach  the  dam  was  not  known  for  Tallahala  Creek  Lake.  For  these 
simulations,  the  depletion  of  D.O.  due  to  travel  time  was  neglected  and 
it  was  assumed  that  the  net  D.O.  contributed  by  inflow  was  equal  to  the 
computed  inflow  D.O.  content. 

9.  Observed  thermal  and  D.O.  profile  data  were  available  for  four 
Mississippi  lakes:  Okatibee  Lake  (1969-71),  Enid  Lake  (1971-72),  Sardis 
Lake  (1968-7I),  and  Grenada  Lake  (1970-72).  Analysis  of  these  profile 
data  yielded  a range  of  values  for  total  depletion  rates,  surface  satura- 
tion percentages,  and  depths  to  which  saturation  extends. 

Dissolved  Oxygen  Simulations 


I 
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10.  Simulations  were  conducted  to  evaluate  D.O.  conditions  within 
and  downstream  of  Tallahala  Creek  Lake.  The  intake  configuration  used 
was  the  type  II,  ns  discussed  previously.  Five  D.O.  conditions  were 


simulated  usitif'  dU'fcrent  depletiojj  I'ates,  surface  saturation  pcrceiit- 


at:’es , and  deptlu;  to  wliicii  saturation  extends.  Tiie  values  used  represent 
tiie  ruripic  of  values  determined  from  tlie  observed  data  on  tiie  Mississipjii 
lakes.  Tiie  input  values  used  arc  sliown  in  Table  Al,  Plate  A?  sliow:: 

D.O.  content  of  flow  entering  t.tie  Intake  structure  for  eacli  of  tiie  five 
study  years  for  all  five  D.O.  conditions.  Plate  A3  sliows  D.O.  profiles 
wltliin  Tallaliala  Ci’eek  Lake  for  eacli  of  ttie  five  study  years  for  all 
five  D.O.  conditions.  Analysis  of  tiie  predicted  D.O.  profiles  (Plate  A3) 
yields  total  depletion  rates  (0.05-0.12)  wlricli  are  witlrin  tiie  ranee  of 
the  observed  data  on  tiie  existine  Mississippi  lakes. 

11.  Tlie  mathematical  model  computes  a value  for  the  D.O.  content 
of  tlie  withdrawal  as  it  enters  tlie  outlet  structure.  It,  is  known  that 
rcaeration  can  occur  as  the  result  of  flow  di scharei ne  under  tlie  reeulat,- 
ine  eate  and  throu/'h  the  conduit,  and  stilline  basin.  Ri'aerat.ion  dat.a 
from  four  noi'thern  Mississippi  impoundment. s (Enid,  Oardis,  Grenada,  and 
Arkabutla  I,akes)  were  readily  available  and  were  used  to  develop  fteneral 
infoniiation  rct'ardinp,  D.O.  uptake.  Plate  Ah  shows  the  percenta/^e  chanr.o 
in  D.O.  sat.uration  t.hrour.h  the  outlet,  works  and  stilliny,  basin  as  a func- 
tion of  the  initial  percentaf^e  of  D.O.  sat.urat.ion  of  the  flow  withdrawn 
from  these  reservoirs.  Plates  A5  and  a6  show  schematic  representations 
of  the  outlet  works  at  the  four  nortiu'rn  Mississippi  imiioundments . 

Plate  AT  shows  a schematic  representation  of  t.he  Tallaliala  Creek  Lake 
outlet  works.  Gince  the  out.let  works  and  stillin;::  basins  at  these 
northern  Mississippi  impounitments  are  (iifferent  from  the  outlet  works 
.and  stillintt  basin  of  the  proposeii  Tallaliala  Creek  liake,  the  data  from 
Plate  Ah  should  not  be  interpreted  as  a prediction  of  reaeration  for 
Tallaliala  Creek  Lake.  Hat. her  these  data  are  intended  to  show  that  the 
D.O.  content  of  reservoir  outflow  can  be  substantially  increased  as  the 
flow  passes  throu/'h  the  out.let.  works  and  stilliiifT  basin,  unless  the  flow 

was  initially  almost  fully  saturated.  Recent  prot,otyi>e  tests  at  Reltx- 

12 

ville  Dam  in  I’ennsylvania  yielded  similar  results,  that  is,  consider- 
able increases  in  release  D.O.  content,  and  additionally  showed  that  the 
predominance  of  reaeratlon  occurred  in  the  outlet  conduit  immediat.ely 
downstream  of  the  w.ater-qual i ty  control  gate.  Calculations  performed  by 


Ah 
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SANT  for  the  Tallahala  Creek  Lake  outlet  works  also  Indicated  that  re- 
aeration of  flow  passing  the  outlet  works  would  be  substantial.  Finally, 
although  many  variables  affect  the  reaeration  potential  of  the  flow,  re- 
aeration is  associated  primarily  wioL  flow  turbulence.  Tizrbulence  is 

characteristic  in  a hydraulic  jump.  A previous  hydraulic  model  investi- 
lU 

gation  conducted  at  WES  indicated  that  a hydraulic  Jump  will  be  pres- 
ent in  the  Tallahala  Creek  Lake  stilling  basin  for  flows  greater  than 
150  cfs,  and  for  lesser  flows  tiic  nydraulic  jump  will  occur  within  the 
outlet  conduit.  Additionally,  the  regulating  gate  will  induce  turbu- 
lence in  the  flow  downstream  of  the  gate,  even  for  very  low  flows. 

Based  on  the  above  considerations,  it  is  expected  that  the  flow  passing 
through  the  Tallahala  Creek  Lake  outlet  works  will  sufficiently  reaerate 
to  p>'oviJe  acceptable  levels  of  D.O.  (5  mg/ 5-  minimum)  in  the  release  im- 
mediately downstream  of  the  structure,  regardless  of  the  D.O.  content  of 
the  flow  entering  the  outlet  works. 
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Table  A1 


Condition 

Depletion 

Rate 

Surface 

Saturation 

percent 

100 

Saturation 

Depth 

1 

0.05 

1°C» 

2 

0.12 

100 

O 

O 

3 

0.12 

100 

5 ft 

1* 

0.12 

80 

5 ft 

5 

0.20 

100 

] ft 

* 1°0  is  the  depth  at  which  the  temperature  is  1°C  less  than  the  sur- 

face temperature. 
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APPENDIX  B:  EFFECTS  OF  A SIMPLE  OPERATION  PLAN 

1.  The  original  simulations  of  Tallahala  Creek  Lake  were  conducted 
by  allowing  the  model  to  select  the  appropriate  ports  during  every  time 
step  in  order  to  best  meet  the  specified  downstream  temperature  objec- 
tive. This  daily  optimal  procedure  is  optimal  in  the  sense  that  the  re- 
lease temperature  for  any. day  is  as  close  as  possible  to  the  objective 
temperature  for  that  day.  However,  a lake  is  a dynamic  system,  and  a 
series  of  optimal  releases  does  not  provide  the  same  results  as  a sched- 
ule of  optimal  releases  based  on  an  entire  year  of  simulation.  For  ex- 
ample, the  release  of  cold  water  to  meet  an  objective  temperature  in  the 
spring  or  summer  might  deplete  the  cold  water  required  in  the  early  fall. 
In  addition  to  the  problems  involved  in  determining  an  operation  sched- 
ule which  will  best  meet  downstream  water-quality  requirements,  the  ac- 
tual operation  of  a lake  according  to  a daily  optimal  schedule  is  not 
necessarily  the  most  practical  mode  of  operation.  Daily  optimal  opera- 
tion may  require  continual  gate  changes  based  on  current  and  expected 
future  temperature  conditions.  Considering  the  greater  requirements  of 
daily  optimal  operation,  a simple  operation  schedule  was  simulated;  and 
the  results  were  compared  with  results  of  simulations  in  which  daily 
optimal  releases  were  used. 

2.  For  these  additional  simulations,  the  intake  configuration 
used  was  type  II  as  discussed  in  the  main  text  of  this  report.  D.O.  con- 
tent was  simixlated  using  a constant  oxygen  depletion  rate  of  0.12  mg/ A/ 
day  and  surface  saturation  of  100  percent  down  to  the  depth  at  which  the 
temperature  was  1°C  less  than  the  surface  temperature.  The  operation 
plan  investigated  consisted  of  the  following  periods. 

a.  1 January-25  June:  release  of  warmest  water  possible;  all 
flow  is  from  the  highest  selective-withdrawal  ports  if 
possible;  excess  is  released  through  floodgate. 

.b.  26  June-27  October:  release  of  coldest  water  possible; 

flow  is  from  the  lowest  selective-withdrawal  ports  except 
for  flood  flows  which  can  be  released  from  the  floodgate. 
This  period  was  included  in  an  attempt  to  delay  the 
buildup  of  water  in  the  bottom  with  low  D.O.  content. 
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2_.  28  October-31  December:  release  of  warm  water  from 

highest  selective-withdrawal  ports.  During  this  period, 
release  flow  will  not  include  low  D.O.  water  from  the  bot- 
tom except  when  flood  releases  are  required. 

3.  Computed  release  temperatures  and  D.O.  concentrations  are  pre- 
sented in  Plates  B1  and  B2.  Comparison  of  these  results  with  Plate  7 in 
the  main  report  and  Plate  A2  (sheet  2)  in  Appendix  A indicate  that  re- 
lease temperature  and  D.O.  content  are  nearly  the  same  for  either  the 
daily  optimal  procedure  or  the  simple  operation  plan.  The  simple  opera- 
tion causes  a large  drop  in  release  temperature  and  D.O.  content  as  a 
resiilt  of  the  change  from  warm-water  release  to  cold-water  release. 
However,  this  transition  could  be  accomplished  gradually,  thereby  reduc- 
ing the  sharp  1-day  temperature  and  D.O.  gradients  for  the  release. 
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APPENDIX  C : NOTATION 


A,B,C,D  Coefficients  of  objective  temperature  curve 
A -13. 83^* 

B 1.721  X 10"^ 

C 1.333 

D 65. hh 

D.O.  Dissolved  oxygen 

D.O..  Inflow  D.O.  content,  mg/ i. 

inflow 

D.O.  Temperature-dependent  saturated  D.O.  content,  mg/ I 

e Natural  logarithmic  base  (2.7183) 

E Mean  daily  equilibrium  temperature,  °F 

2 

H Net  rate  of  surface  heat  transfer,  Btu/ft  /day 

2 

Rate  of  heat  absorbed  in  layer  (i),  Btu/ft  /day 

H Rate  of  heat  transfer  into  or  out  of  surface  layer, 

^ Btu/ft^/day 

2 

K Coefficient  of  surface  heat  exchange,  Btu/ft  /day°F 
Q Mean  daily  streamflow,  cfs 

2 

5 Rate  of  total  incoming  shortwave  radiation,  Btu/ft  /day 

t Julian  day 

T Stream  temperature,  °C 

Depth  below  surface,  ft 
a,B  Regression  coefficients 

a 19.^75 

Mixing  coefficient  at  surface 
Mixing  coefficient  at  bottom 

6 Percentage  of  incoming  shortwave  radiation  absorbed  in  the 
surface  layer 

6^  -0.0020 

0.13595 

0.123*t 

O.U095 

9 Mean  daily  stream  temperature,  °F 

X Absorption  coefficient,  ft”^ 

Cl 
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